:j0_k_norm => (log10 => L"Log Normalized Current Density ($J_A$)")
In [2]:
Code
begin j_events_low_fit =load() j_events_tau_20_fit =load(tau =20) j_events_high_fit =load(ts =0.12) dir ="../data/05_reporting" w_events =load("$dir/events.Wind.fit.ts_0.09s_tau_60s.arrow")# add a label column to the dataframes j_events_low_fit[!, :label] .="1 Hz (fit)" j_events_high_fit[!, :label] .="8 Hz (fit)" j_events_tau_20_fit[!, :label] .="1 Hz, 20s (fit)"# filter high time resolution events j_events_high_fit =@chain j_events_high_fit beginfilter(:len =>>(240), _)end# combine the dataframes j_events =reduce(vcat, [j_events_low_fit, j_events_high_fit, j_events_tau_20_fit])println("Number of events: ", size(j_events, 1))end
┌ Warning: automatically converting Arrow.Timestamp with precision = NANOSECOND to `Dates.DateTime` which only supports millisecond precision; conversion may be lossy; to avoid converting, pass `Arrow.Table(source; convert=false)
└ @ Arrow /Users/zijin/.julia/packages/Arrow/Y6R1E/src/eltypes.jl:273
Number of events: 61879
In [41]:
Code
functionload_tau(tau) df =load(tau=tau) df.label .="$tau s"println("Number of events: ", size(df, 1)) dfendj_events_taus =60:-10:20.|> load_tau |> x ->reduce(vcat, x)
Number of events: 20875
Number of events: 22721
Number of events: 24788
Number of events: 28132
Number of events: 32944
# plot the number of events with respect to tau in different radial distancesspecs =data(j_events_taus) *mapping(:tau, row=:r) *histogram()draw(specs, facet=(; linkyaxes=:none))
In [47]:
Code
# get unique events defined by the starting time and ending timecols = [:"t.d_start", :"t.d_end"]println("Number of events: ", size(j_events_taus, 1))j_events_taus_u =unique(j_events_taus, cols)println("Number of unique events: ", size(j_events_taus_u, 1))
Number of events: 129460
Number of unique events: 88679
Check the discontinuities properties with radial distance
# groupby r and describe the data for each group # j_events |> @groupby(_.r) |> @map({r=key(_), j0_k=describe(_.j0_k), L_k=describe(_.L_k)})@chain j_events begingroupby(:r)combine(:plasma_density => mean, :ion_inertial_length => mean, :b_mag => mean) end